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Systems with topologically protected ground-state degeneracies are currently of great interest due to their 
potential applications in quantum computing. In practise this degeneracy is never exact, and the magnitude of 
the ground-state degeneracy splitting imposes constraints on the timescales over which information is topologi¬ 
cally protected. In this Letter we use an instanton approach to evaluate the splitting of topological ground-state 
degeneracy in quasi-ID systems with parafermion zero modes, in the specific case where parafermions are real¬ 
ized by inducing a superconducting gap in pairs of fractional quantum Hall (FQH) edges. We show that, like ID 
topological superconducting wires, this splitting has an oscillatory dependence on the chemical potential, which 
arises from an intrinsic Berry phase that produces interference between distinct instanton tunneling events. 
These Berry phases can be mapped to chiral phases in a (dual) quantum clock model using a Fradkin-Kadanoff 
transformation. Comparing our low-energy spectrum to that of phenomenological parafermion models allows 
us to evaluate the real and imaginary parts of the hopping integral between adjacent parafermionic zero modes 
as functions of the chemical potential. 
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The possibility of creating bound states with exotic non- 
Abelian statistics at the ends of quasi-ID systems is an ex¬ 
citing prospect that has been developed both theoretically 1-16 
and experimentally 17-2 '' over the past few years. These bound 
states are known to generate a topological ground-state de¬ 
generacy. In arrays of ID systems, states in the resulting 
low-energy Hilbert space can be entangled and manipulated 
by braiding processes 26 , which are (at least in principle) 
robust to noise. Such systems have drawn significant in¬ 
terest due to their potential to realize topological quantum 
computation 27- ' 1 . 

Parafermion bound states 6-16 , which have k ( k > 2) topo¬ 
logically degenerate ground states, are particularly promis¬ 
ing for such quantum computing applications: In compari¬ 
son to Majorana bound states (k = 2), parafermion bound 
states allow for a denser (albeit non-universal) set of com¬ 
putational gates, and are believed to be intrinsically more ro¬ 
bust to environmental noise 8-10,32 . Parafermions are signifi¬ 
cantly more challenging to produce than their Majorana coun¬ 
terparts: Most proposals entail generating them as defects in 
2D Abelian FQH states, in ways that have yet to be carried 
out experimentally. However, their non-Abelian statistics are 
also more complex than for Majoranas, which makes realizing 
them a particularly exciting prospect. 

For bound states confined to the endpoints of a system of fi¬ 
nite length L, the topological ground-state degeneracy is split 
by an amount A E ~ e~ L ^, where f is proportional to the 
correlation length in the bulk of the system, causing superpo¬ 
sitions of ground states to decohere over time. Though the¬ 
oretically it is possible to make this splitting as small as re¬ 
quired by making L large, there are definite advantages when 
this splitting can be made small even for modest-length sys¬ 
tems. For ID topological superconducting wires 33,34 or spin- 
Hall based superconductors 2,5 , this can be achieved by small 
adjustments in the appropriate chemical potential, since in 
addition to the exponential fall-off in L the splitting has an 
oscillatory dependence on the Fermi momentum via A E ~ 
e~ L cos(fcpL). The observed oscillations of the splitting of 


the zero-bias conductance peaks 20 can be viewed as evidence 
of Majorana bound states 3 ' - ’ 7 . 

Despite the publicity it has had in Majorana systems 38 ^ 11 , 
little attention has been given to the ground-state splitting in 
parafermion systems. Specifically, one might wonder whether 
these exhibit an analogue of the oscillatory cos(& fL) term. 
In this Letter, we calculate the splitting of the topological 
ground-state degeneracy in parafermion platforms obtained 
by inducing superconductivity or ferromagnetism at certain 
types of edges in FQH states 8-13 . Following Ref. 34 in the Ma¬ 
jorana case, we perform this calculation using a bosonised de¬ 
scription of the strongly interacting ID system, in which the 
splitting of the ground-state degeneracy is obtained by an in¬ 
stanton calculation in the resulting sine-Gordon model. Inter¬ 
estingly, as in the Majorana case we do find oscillations in the 
splitting as functions of chemical potential or applied mag¬ 
netic field. These intriguing oscillations result from a Berry 
phase term in the generic sine-Gordon action. We also use our 
calculation to deduce the magnitude and phase of the hopping 
coefficients that arise most naturally in parafermion chains, 
such as those studied by Refs. 6 - 12,42 ^ 15 . 

Model of parafermion zero modes .— Several groups 8-13 
have suggested that parafermion zero modes can be generated 
in systems with counterpropagating chiral edges separating 
two FQH regions with opposite ^-factors. The edge of interest 
consists of one right-moving and one left-moving mode with 
opposite spin polarizations. Two types of electron tunneling 
processes can open a gap at this edge: Inducing superconduc¬ 
tivity (SC) generates a Cooper-pairing A(i/r L ff'p - + H.c.), 
while spin backscattering B(ff L fipRA + H.c.) can be induced 
by tunnel-coupling the edge to a ferromagnet (FM). Any in¬ 
terface between these different types of induced gaps will host 
parafermion bound states (Fig. 3). 

To describe this system, it is convenient to bosonise the two 
edge modes, representing the right- (left-) moving electrons 
as f>' R / L ~ ^ 2 ! W ^ ere Pr/l are chiral bosonic 
fields, and 1/n is the filling fraction of the corresponding FQH 
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FIG. 1: (color online). Schematic spatial profiles of the proximity- 
induced gaps A(x) and B(x) for the FM-SC-FM setup. 


bulk regions. Here <; is related to the inverse energy cutoff 
(^cutoff) of the bosonised theory via g ~ vh/E c uto ff. where v 
is the velocity of the edge modes. 

The two backscattering terms are most simply expressed in 
the basis </> = + <p L 4 ) and 0 = ~ W;)- 

These non-chiral fields are related to the charge density pc 
and spin density ps via pc = ^d x 9 and pg = \d x (j). In this 
bosonised basis, the two backscattering terms take the form 
A((/>]) i+H.c.) ~ A sin(2 n<f>) and B(^ L +H.c.) ~ 

Bsm(2n6), and the ID parafermion system is described by 
the following Euclidean sine-Gordon action 1J ' 16 !: : 

Se = [d,Tdx{ih—d x 9(x , t)g> t </(x, r) — ^ d x 9{x, r) 

J n 7T 

+ ~^—{9 x 9(x, r )) 2 + ^^[sin( 2 n 0 (x, r)) + l] 

Z7T 7TO£ 

+ ^-(9 x ( I ) {x,t )) 2 + —— [sin( 2 n 0 (x,T)) + l]}. ( 1 ) 

Here // represents the chemical potential and B. A are energy 
gaps induced by spin- and charge-backscattering processes, 
respectively. The commutation relation [</(x), 0(x')] = 

i^O(x — x') dictates that only one of the two possible gap¬ 
ping terms can have a non-vanishing expectation value at a 
given spatial position. However, if a region where A ~ hv /y 
can be sandwiched between two regions where B ~ hv /<,, 
parafermion bound states arise at the interfaces between them. 
In the bosonised picture, the resulting topological ground- 
state degeneracy is manifest in the 2 n values of <f> for which 
the sine term is minimized 7-12 . 

In the following, we consider the FM-SC-FM heterostruc¬ 
ture shown in Fig. 3, on which we take B = 0 (A, p, = 0) 
in the SC (FM) region \x\ < L/2 (|x| > i/2). In the FM 
regions \x\ > i/2, the field 9 is therefore pinned to one of 
the potential minima, and <j> is strongly fluctuating as required 
by the commutation relations. As we show in A, under these 
conditions the FM regions do not contribute to the ground- 
state energy splitting, and after integrating out 9 we obtain the 
following effective action for the SC region: 

TV L 

S 4 , = [ drf dx{^(d T (j)(x,T )) 2 + ^-(d x tj>(x, t )) 2 

H-[sin (2 n<j>(x, r)) + 1 ] + i^—^d T (j){x, r)}. ( 2 ) 

7r ng 7 tv 

In the ground states of this reduced system the </> field is also 
approximately pinned at one of the 2 n inequivalent local min¬ 
ima <p m i„ of the sine potential; we will take A to be suffi¬ 
ciently large that these low-energy states are well separated 


from the rest of the spectrum. The fluctuation-induced split¬ 
ting between the 2 n otherwise degenerate ground states is then 
determined by the amplitude for tunneling between adjacent 
local minima. 

The last term in Eq. (2) plays the role of a topological Berry 
phase term ,S‘b- p , contributing to the net action only for field 
configurations which start and end at different values of 9 (i.e. 
only for instantons). It introduces oscillations in the splitting 
of the ground-state degeneracy as the chemical potential // is 
varied. To the best of our knowledge, 5 b- p was not included in 
previous studies of instantons in the bosonised periodic sine- 
Gordon model, which considered the case // = 0 34,48,49 . 

Instanton calculation of level splitting .— For the sine- 
Gordon model described by Eq. (2), the minima of the poten¬ 
tial term are at <j> min = - s ^ where j = 0, 1 , ..., 2 n - 1 , 
with the j- th ground state denoted by | j). The classical 
soliton solution interpolating between vacua at j — — ■£■ and 
(j ± 1 )))—// has the form 50-52 : 

+j-±- arctan[e“ (T_ro) ], (3) 
4n n n 

where u> = 2y/Av/(hg) >0. To a good approximation, we 
may neglect spatial variations in the instanton solution due to 
boundary effects (see Appendix A). 

Following Refs. 49,53-56 , the amplitudes for the one- 
instanton and one-anti-instanton processes are: 

{j\e~ H * Tr/h \j T l)o.i. 

(4 » 

Here So + (—)ih , y is the effective action of an instanton (anti- 
instanton) tunneling event (with small quantum corrections 
omitted, as discussed in Appendix C), where the imaginary 
contribution 7 stems from Sb-p, and hoj /2 stands for the zero- 
point energy of the ID harmonic oscillator. In terms of the 
parameters in Eq. (2), we have So = , 7 = . As antici¬ 

pated, the amplitude for tunneling vanishes exponentially with 
the wire length L provided the bulk is gapped (i.e. A f 0). 
More significantly, we now perceive the importance of the 
Berry phase term, which contributes a different net phase to 
the amplitude for instanton and anti-instanton processes. 

The prefactor K, = ^vuiSo/(2nhL) is the Fredholm de¬ 
terminant describing Gaussian fluctuations about the saddle- 
point solution (3). This determinant is sensitive to spatial fluc¬ 
tuations in the SC region, and its scaling with L is sensitive 
to the choice of boundary conditions (BCs). As shown in 
Appendix C, for Neumann-type BCs appropriate to the setup 
shown in Fig. 3, this prefactor is independent of L, in agree¬ 
ment with existing work on Majorana bound states 33,35 57 . For 
the energy splitting in periodic BCs relative to antiperiodic 
ones, on the other hand, the prefactor contains an extra factor 
of Fig. 2 compares this prediction to numerical val¬ 

ues for the energy splitting in the Majorana nanowire for both 
BCs. 

Armed with one-instanton solutions, the total instan¬ 
ton contribution is obtained using a dilute instanton gas 
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where a denotes the parafermion bound states, and F is a 
part of the topological data. For the case at hand the possi¬ 
ble choices of F are described in Ref. 59 ; taking [F““] ao = 
e i(7T/n)[(a-g*) mod 2n\ and = -h/Ce- s °/ h -^ 6 a p recovers 
the form (7) for the energies of these ground states. 

It is instructive to check Eq. (7) for the case of Majorana 
fermions (n = 1), where the splitting can be calculated di¬ 
rectly from a quadratic fermion Hamiltonian 5 ’ 33,38 - 39 . The rel¬ 
evant calculation in the quantum Hall systems described here 
can be carried out as for the nanowire case 33,35 - see Ref. 97 
and Appendix B. The result is 

A B _a l ( /iL\ 

E “ = - C ATB e ~ hVCOS V q -J^)' 9 = °’ 1 ’ (9) 


FIG. 2: (color online). Numerical results for the oscillatory zero- 
energy splitting as a function of L in the ID nanowire models hosting 
Majorana bound states 1 ' 3 ' 4 , for (a) open and domain-wall BCs, and 
(b) periodic versus antiperiodic BCs (PBC-APBC). (Recall that for 
a Majorana SC ring, the ground state with PBC (APBC) is a parity- 
even (parity-odd) state.) The term linear in L has been subtracted off 
in both cases, leaving a result independent of L in (a), and depending 
logarithmically on L in (b). 


approximation 34 , which gives: 

(j + \e~ H ^ /h \j-) =Af [ 
Jo 


d( 

2 t/ 




x exp 


-~T t + 2K.T T e~ So/h cos(C - 7 ) 


(5) 


The cosine term can be viewed as arising from interference 
between instanton and anti-instanton trajectories, for which 
S'b-p has the opposite phase. 

In the presence of instanton tunneling events, the eigen¬ 
states of H, : , are therefore Bloch-wave-like states of the form: 

|C) °C (^) 5 e_i ° U>- Imposing (j + ±2n\e~ H * T ^ n \j- ± 
2 n) = {jj r \e~ H,l,TT ^ h \jJ) to account for the fact that only 2 n 
of these minima are physically distinct forces ( to take on the 
2 n discrete values -q^, with q^ = 0,1,..., 2n — 1, which 
gives the 2 n states: 


M oc 



2n—1 


E e ~ ilw \ o- 

3=0 


( 6 ) 


From Eq. (5), the energies of these states are, up to an overall 
constant. 


E (q<») 


2 hu! 

IT 



7r pL \ 

n hnv J ’ 


(7) 


where g 0 is the generalised charge parity conjugate to the dis¬ 
crete, compact variable © mm . This form agrees with the gen¬ 
eral result of Ref. 78 for the energies of anyon pairs: 


2n-l 

E(q<p)= E (ra[F“ aa ]^+H.c.), 

a —0 


where C is a constant of order unity. The coefficient of the 
decaying exponential differs from Eq. (7), since for n = 1, 
Sq/H = (AT) ^ Hv / (As). However, in the instanton cal¬ 
culation of the exponential term we neglect all modes above 
the gap set by A; hence the cutoff energy hv/s should not be 
much larger than A. Further, the tunneling process requires 
the presence of virtual fluctuations up to an energy of ap¬ 
proximately A, so the cutoff energy should also not be much 
smaller than A. The factor - sJhv/(As), which parameter¬ 
izes the difference in the exponential decay lengths from the 
two calculations, is therefore a constant of order 1. Parallel 
reasoning applies to the comparison of the prefactors. Re¬ 
markably, the argument of the cosine term agrees exactly with 
our instanton calculation, suggesting that this oscillatory de¬ 
pendence on /j is insensitive to the BCs and to the various 
approximations being made. 

We note that the preceding analysis also applies to the case 
B > 0, A = 0 by taking <j> ^ 0 and replacing the chemical 
potential term with a magnetic field term of the form —yd x cf>. 
Therefore the oscillatory dependence of the ground-state en¬ 
ergy splitting on chemical potential or magnetic field is a rel¬ 
atively ubiquitous feature of parafermion zero modes. 

Hopping in parafermion chains. —One interesting applica¬ 
tion of our calculation is that it allows us to infer the phase 
of intrawire parafermion hopping terms. This is of particu¬ 
lar interest as chains of coupled parafermions can be used to 
generate even more exotic topological phases 12,44 . 

In the setup we consider, the parafermion bound states can 
be described by operators a^\a^ which annihilate (create) 
parafermion zero modes at the left and right endpoints of the 
SC region, respectively, and satisfy the relations 

a l/R = 1) a L/R = a ( l}r and a L&R = a R a L e l », (10) 

which are sufficient to ensure that these bound states have 
“parafermionic” non-Abelian statistics 7-1 1 . 

In terms of the bosonised fields 9 and we have a\a.R = 
e l ™( ?-1 / 2 ), where q = ^(9(j) — 9(— ^)) is the total charge in 
the SC segment modulo 2 s 1 . From the commutation relation 
between 0 and 9 , it follows that 



( 8 ) 


































4 


which carrying out precisely the tunneling processes whose 
matrix elements we have just evaluated. 

The low-energy Hamiltonian describing the parafermion 
tunneling between the two endpoints is therefore 

H a = ta' L a R + t*a R a L - 


Its eigenstates are labeled by an integer q = 0,1,2,..., 2n— 1, 
and satisfy 8-11 

ot\oL R \q) = —e z " ( «“5)| q), a ] R a L \q) = -e - *" (9_ a) \q). 

( 11 ) 

Note that Eqs. ( 6 ) and (11) together also fix the phase associ¬ 
ated with the action of the parafermion hopping term on the 
eigenstates of (j>\ a\a R \j) = j + 1 ), a* R a L \j) = 

e + ’' 2 ™ |j — 1). The corresponding energies—which are pre¬ 
cisely the energies that we have just obtained with our instan- 
ton calculation—are: 


E(q) = -2^/t| +f|cos 


( 12 ) 


where we have defined tan$ = fcj/tsR. 

Comparing Eqs. (7) and (12) allows us to constrain £sr and 
£q . For Majorana fermions (i.e. n = 1), there is an additional 
constraint: Since = a R / R , the two terms in Eq. (11) are 
not independent. This forces t sr = 0 (i.e. d = ±7r/2), and 
fcj = ±hK.e~ s °^ n cos [[iL/ ( hnv )]. For n > 1 there is no such 
a constraint; however, in these cases matching the eigenvalues 
of both H and the operator a\^a R = e l ( e (^) _0 ( _ ' J ) _,r /( 2n )) 
fixes d, such that 

tsR = ±h/Ce~ s °/ h cos[7r/(2n) — /.iL/(hnv )], (13) 

fcj = ±h/Ce~ s °/ h sin[7r/(2n) — fiL/(hnv)\. (14) 

Using the analogous approach for an SC-FM-SC system (with 
[j, => h, A => B) gives the analogous conclusion. 

It is worth stressing that even at vanishing /i (or h), for n > 
1 the hopping parameter t is complex with d = 7r/(2 n). This 


suggests that the proposal for universal quantum computing 
by manufacturing Fibonacci anyons in coupled parafermion 
chains l: - 42 - 44 - 45 is better achieved in systems without finite 
chemical potential or magnetic field. More specifically, the //- 
dependent contribution to d in Eqs. (13) and (14) corresponds 
to a chiral phase 60-63 in the quantum clock model. To be con¬ 
crete, a system of 2 N tunnel-coupled parafermion zero modes 
is dual (via the Fradkin-Kadanoff mapping 6 - 46 64 ) to an 7V-site 
chiral quantum clock chain, where with appropriate conven¬ 
tions hopping across a SC (FM) region maps to the trans¬ 
verse field (ferromagnetic clock) coupling. Under duality, the 
phases of the parafermion hopping terms t map to a chiral 
phase of ±e* for the ferromagnetic clock coupling, as well 
as a chiral phase of ±e 1 '^ for the transverse field term. No¬ 
tice that the oscillatory dependence of the ground-state degen¬ 
eracy splitting on the chiral phases and the system’s size has 
been observed numerically for these chiral clock systems 45 
recently. 

In summary, our non-perturbative calculation shows that it 
is possible, in principle, to tune the magnitude of the ground- 
state splitting in parafermion systems (as well as the phase 
of the parafermion hopping parameter, for n > 1 ) by means 
of a chemical potential or an external magnetic field due to 
interference between distinct instanton trajectories resulting 
from a topological term in the effective action. Because the 
period of the resulting oscillations is given by giL/(hnv) (or 
hL/(hnv)), this splitting can be fine-tuned with relatively 
small changes in //. (or h). As for Majoranas, we anticipate 
that this fact will be both of practical use to achieve quantum- 
coherent systems, and a potential signature of the existence 
of parafermions in these systems. Finally, our results might 
also be applicable to the spin-unpolarized v = 2/3 FQH het¬ 
erostructures proposed by Refs. 1213 . 
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Appendix A: The reduced effective action for instanton configurations 

Here we justify the effective action S 0 which is our starting point for studying instanton effects in the main text. Specifically, 
we will show that, for a short SC region sandwiched between two long FM regions, the effect of instantons is limited to the 
SC region. We will use this to derive an effective action in the SC region alone, with appropriate boundary conditions at the 
boundaries of this region. 

To address the details of this problem, we must be slightly more precise in what we mean by FM and SC regions. We 
consider the FM-SC-FM heterostructure where the edge is gapped by FMs when x < — and x > , while in the interval 

- < x < , the edge is gapped by SC. To simplify the calculation, we approximate these gaps as step functions (see 

Fig. 3). The commutation relations between 0 and <f> make it impossible to induce both kinds of gaps at the same spatial position; 
therefore there is necessarily a small domain wall (DW) region between FM and SC regions, whose width we take to be 5 <C L. 

Our derivation proceeds in three steps. First, we will argue that, provided the FM gap is large, the instanton configurations of 
interest have essentially no effect on physics in the FM region. This allows us to restrict our attention to the effective action in 
the SC and DW regions, with appropriate boundary conditions at each domain wall. Second, we integrate out 6 in the SC and 






5 


DW regions to obtain an effective action for the field </>. Finally, we argue that the specified boundary conditions mean that the 
contribution of the DW region to the tunneling amplitude is also essentially independent of the instanton configuration, allowing 
us to study the ground-state splitting using an effective action for the SC region alone. 



L+S L-S L-S L + S X 

2 2 2 2 


FIG. 3: (color online). The spatial profiles of the proximity-induced gaps A (a;) and B(x) for the FM-SC-FM setup. The two gapped regions 
are necessarily separated by a domain wall region of width 5 that accommodates a parafermionic zero mode denoted by the green ellipse. 


The FM region. First, we will argue that the impact of the FM region on the ground-state splitting can be captured simply 
by choosing appropriate boundary conditions for the fields 9, tf> at x = ±^y^. Essentially, this is because the instanton profile 
vanishes for |.t| > , and hence the path integral on these regions is indifferent to instanton effects. 

To see why this must be so, we begin by considering the limit yy 00 • In this limit 9 is exactly frozen in each FM region to 
one of its 2 n minimal values, and its spatial and temporal derivatives in this region vanish. Since [— d x 9(x), </>(x)\ = i, pinning 
d x 9 = 0 in this region forces <j> to be evenly distributed on the interval [0, 2i r), irrespective of the behaviour of <f> for \x\ < . 

Hence in the limit that 9 and its derivatives are perfectly pinned in the FM regions, their contribution to the tunneling amplitude 
is necessarily independent of the instanton configuration in the SC and DW regions. Hence in this limit, we can replace the 
FM region with free boundary condition for cj>(±(L + 5)/ 2), and a sum over all possible pinned values of 9 at each of these 
boundaries. (This sum ensures that the system can access all possible values of q.) 

In other words, for -¥ 00 the tunneling amplitude for the entire system can be expressed: 

X= ^ f d(j)id(j)2ZFu[9L, 4 >i]T[0l, 9r, <j>u <fa\ZFM.[0R, <fa\, (Al) 

Sr,6 l J 

where the summation is over all values of 9r,9r £ -{0,1,..., 2n — 1}. Here we have defined </>(—(L + 6)/ 2) = 
4>i, (f> ((L + 5)/ 2) = (f> 2 . Note that we must integrate over all values of <j>i, </> 2 , since d x 9 is pinned at this boundary. 

From the form of the action in the FM regions (which depends only on derivatives of <j>, with no potential term), it is evident 
that Zpm will be independent of the particular values of and also of the choice of 9l,b. within the set of minima of the 
potential. Therefore any impact from instanton configurations must be contained in the tunneling amplitude T[- ■ ■ ], which 
describes the SC and DW regions. 

Hence provided B is sufficiently large that fluctuations in 9 in the FM region can be neglected, the FM region cannot contribute 
to the instanton transition amplitude. Including small fluctuations, to account for the fact that d x 9 is not exactly pinned in the 
FM region, we would find that the boundary values of (j> need not be perfectly uniformly distributed, making some “leaking” 
of the instanton configuration into the FM region possible. However, provided the fluctuations in 9 are small (i.e. the FM gap 
is large), the distribution of (!> must be wide. The finite correlation length in this region also ensures that any correlation with 
the instanton solution occurs over a finite (and small) spatial extent. Therefore though the impact of the instanton configuration 
on the FM region in this case is not strictly 0, we will neglect it relative to the remaining terms, an approximation that is valid 
provided the FM gap remains large. Including such terms could result in corrections of order 1 of the magnitude of the instanton 
transition amplitude, but will not contribute to the Berry phase term in the action provided /r = 0 in the FM regions. 

Effective action in the SC and DW regions. Based on the previous discussion, to calculate the ground-state splitting we 
may focus on the region \x\ < y^A, where the potential term for 9 vanishes. The boundary conditions on this region are 
that (j) is free at both spatial boundaries, and that 9 (— = 0l, 9 (yy) = 9r, where the transition amplitude contains a 

sum over the possible values ^9rr = 0,1,. . . , 2n — 1. Recall that the total charge in the SC segment is fixed by q^jn = 
n (0 (y^) — 9 (y yy)) = 0,1/n,..., (2 n — l)/n (mod 2); hence this sum ensures that all possible charges are included in 
the path integral. 
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Within this region, the potential term for 9 vanishes, and we may complete the square in the path integral to obtain: 

2 


S= dr 


l\x\< 


hnv 
dx < —— 

l+s 27 r 


d x 9{x,r) - + -d r (j>(x,T) 

hnv v 


hn /rx f/ hnv A(x) r . 

+ x— (a T 0 (x,r))" + — (d x <j>(x,r)) H-- sin( 2 n^(x,r)) + 1 ] 

Z7n/ Z7r 7rn^ 

, .•(‘Wa A (m (^)) 2 

+ *-Or</>(a;, t) - —-— 

7 T 7 / An hnv 


(A2) 


The second line in Eq. (A2) is the standard sine-Gordon action 5 s _g, and the first term in the third line is the newly-derived Berry 
phase term Sb- p . Shifting 9 and performing the Gaussian integral leads to the effective action given in the main text: 


SU = 


/ dr dx{ 

J\x\<±P- l 2 ™ 


hn ,. NX o hnv A (a?) r . /rx / xx 

(5 r 0(x,r)) + (d x (/)(x,T)) + [sm(2 n<j>(x,r)) + 1] 


27r 


7rn^ 


i^-d T <t){x,T) 

7TV 


(A3) 


where we have dropped the constant term. 

In the main text, we ignore the domain wall region, and treat the classical instanton field as position-independent. Relaxing 
this approximation will result in corrections of order 6 /L to the classical instanton action. In particular, the contribution of time 
derivatives of <f> in this region is clearly small in S/L. Further, after averaging over all boundary values at x = the set 
of differences of <j> across the spatial extent of the DW region is independent of the pinned value of © in the SC region (and its 
time derivative) - i.e. the distribution of d x cf> in the DW region is insensitive to the instanton configuration, and hence does not 
contribute to the instanton transition amplitude. 

Boundary terms from integrating by parts. To obtain the action Se given in the main text, we have integrated a term of the 
form d T 9d x 4> by parts. Here we verify that the resulting boundary terms do not contribute to the effective action of the instanton. 
To show this, we use periodic boundary conditions in time for the fluctuating components of our fields. 

|7V /2 1 
I-TV/2 


The boundary terms in question are: J dx ( iO(x,T)d x <j>(x,T )) 


— f dr (i9(x,T)d T (f>(x,T)) 


We begin with the 


first term. In the SC region, we may separate (j) = <fi so i + 77 , where t](x, T t /2) = 77 ( 2 , — T t /2) and d x 4>soi = 0. In the FM region 
<fi = 77 , and since there are no instantons in 9 the configurations at ±T t /2 are identical. In the DW region the field does vary 
spatially, but this variation is not due to 4> so i, but rather due to the fluctuations (in other words, the distribution of d x <p in the DW 
is time-independent, as emphasized above). Therefore the first term vanishes. 

For the second term, (which we compute for the whole system, not only the region relevant to the instanton action), we obtain: 


dr 


(i0(x,T)d r <f>(x,T)) r ~ = 0 R / dr [id T (j)(oo,T)\ - 0 L dr [id T <j)(-oo,T)] 


= 0 r (<K 00 , T r /2) - <(>( 00 , —TV/2)) - 0 L (0(-oo, TV/2) - ^(- 00 , -TV/2)) = 0, 


where we have used the fact that since each FM region is semi-infinite, there are no instantons in 9. The last equality follows 
from the fact that in the FM region </> consists only of fluctuations, which obey periodic boundary conditions in (imaginary) time. 

Thus both of these boundary terms do indeed vanish. 

Shifting 9. Here we present the details of the shift in 9 , to confirm that the resulting integral is independent of the configuration 
of cj). For a one-instanton contribution, we have </>(x,r) = (f>soi{x,T) + 77 ( 2 , 2 ), where the fluctuation field 77 obeys 77 ( 2,2 = 
±TV/2) = 0, and (f> so i is the classical instanton solution, which we argued above vanishes at 2 = ■ We must now evaluate 


V9{ 2 , 2 ) exp 




dx[d x 0(x , 2 ) 


hnv 


+ -d T r)(x,T) + -d T (j) so i(x,T)Y 

V V 


(A4) 


Since the region in question is finite, one might worry that such a shift will result in different boundary conditions for 9 and 9, 
either in space or in time, leading to a different set of eigenvalues even for the same differential operator. However, the boundary 
conditions are not affected by (f> so i or //. First, both of these quantities vanish at the spatial boundary. Second, in the limit of 
interest T t —» 00 , d T (j> so i vanishes at the time-like boundaries, and 77 is time independent. This leaves the fluctuation field 77 , 
which obeys free rather than fixed BCs at the spatial boundaries (but the same BCs in time). This change in spatial BCs between 
9 and 9 does not, however, result in any additional dependence on 77 in the functional integral. 

To perform the desired shift, we Fourier transform all fields by using 0 ( 2 , 2 ) = 5 /-^ YZk iu> h{k,iLo)e lkx e~ lulT (and sim¬ 
ilarly for other fields), where to are real, discrete frequencies, then shift the Fourier components 0(k, no) to absorb the 
chemical potential and the soliton terms into the Gaussian integral of 9. This gives: 0(k, iu>) = 9{k , iui) — ^ ^“ |Q + 
(77 (fc, iu j) + (j> so i(k, iui)), where i s the Kronecker (5-function. The resulting Gaussian integral over 9 is independent of 
the specific form of both </> so i and 77 , leading to the effective action (A3). 
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Appendix B: Splitting of ground-state energy in spin-Hall edge: A wavefunction-overlap calculation 

As an important check, let us discuss the Majorana zero modes at the edge of a 2D quantum spin-Hall system which is in 
close proximity to an s-wave superconductor and two ferromagnetic insulators. We assume that the right movers carry spin up 
and the left movers carry spin down. The Bogoliubov-de Gennes (BdG) Hamiltonian for the edge states is as follows, 

U(x) = -iv F T z a z d x - h(x)t z - A(a ;)t x + B(x)o x . (Bl) 

The Pauli matrices r and a act on the particle-hole and spin spaces, respectively. To focus our discussion, we consider the 
FM-SC-FM setup shown in Fig. 3. Neglecting the width of the domain wall region, this gives the following spatial profiles for 
the chemical potential and the proximity-induced gapping terms: fi{x) = /iQ(x+ ^)@(—x + ^), A(x) = A0(a:+ ^)0 (—x + 
|r), B(x) = B( 1 — 0(x + ^)0(— x+ j)). Here the potential strengths A, B |/i| ^ 0 and 0(at) is the Heaviside step function. 
We also assume that the chemical potential term can be generated by locally gating the edge. 

With the linear dispersion relation in Eq. (Bl), analytical solutions for the Majorana bound states can easily be obtained 
in the limit that the two bound states are infinitely far apart. Specifically, if we divide T~L(x) into the left and right parts for 
x £ (— 00 , 00 ), then H(x) = 'Hl{x) + Ur(x) = 'Hr{x) + Ul(x), where 

/ -iv F d x - Hl(x) B l {x) -A l (x) 0 

Bl(x) iv F d x ~n L (x) 0 -A l (x) 

-A l (x) 0 iv F d x + n L {x) B l (x) 

V 0 -A l (x) B l (x) -iv F d x + h l (x) 


n L {x) = 


(B2) 


U R {x) = 


/ - iht(x ) B r {x) ~A r (x) 0 ^ 

Br{x) —Jxr(x) 0 -Ar(x) 

-A R (x) 0 Mr(x) Br(x) 

V 0 -Ar(x) Br(x) p, R (x) ) 


(B3) 


with n L { x) = fj,0(x + §), A l (x) = A0(r + §), B L (x) = BQ(-x - f), p, R (x) = -fiQ(x - §), A R (x) = 
-A0(i — ^), Br{x) = BQ{x — j). Similar forms can be constructed for Hr(x) and Ul{x) by noticing that Hl{x) = 

KHr(-x)K , U R {x) =U l (-x), Hl(x) = hr{-x), P-l(x) = p, R (-x), A L (x) = Ar(-x), A L (x) = Ar(-x), B l (x) = 
Br(—x). Typically, TLl,r(x) and Ur F {x) are spatially well-separated from each other, which allows to firstly approximately 
focus on the Majorana zero mode from only the part of T-Ll,r{x), then adding the treatment of Ur.l(x) as the perturbation. 

Left Majorana zero mode. We are now in a position to find the analytical wavefunction for the normalizable Majorana bound 
state in the left domain wall x = — 4 by solving the following reduced ID Dirac equation with the proper boundary conditions. 


/ 

-iv F d x - hl(x) 

B l {x) 

-A l (x) 

0 \ 

/ u(x) \ 

B l {x) 

iv F d x - li l {x) 

0 

~A l {x) j 

w(x) 

V 

—A l (x) 

0 

iv F d x + hl{x) 

Bl(x) 1 

w*{x ) 

0 

-A l (x) 

B l (x) 

-iv F d x + hl{x) J 

\ -u*(x) j 


= 0 . 


(B4) 


The self-conjugate Majorana form of the spinor here is fixed by particle-hole symmetry, defined in real space by the operator 
S = r y (g> o y K, where t v ® a y is the charge conjugation matrix and K denotes complex conjugation. 5 anticommutes with 
H(x), and hence constrains the form of the zero-energy solutions. 

T ( 'U'(x') \ 

In the region x > —jj, solving Eq. (B4) yields the 2-component zero-mode eigenvector, Xo( x ) = ( u; * J = 

X_ 1 . ^ \ ^ exp(i^{x + 2 ) — ^(x+ 2 )^, where A > 0, and x~^ is a constant. Similarly, we can solve for the wave- 
' the It 

u ( x ) \ ( 1 ) ex p^Ji_( x _f_ L)^j, with B > 0, and C _l is a constant. The complex constants x _l and 


Co (a) = 


w(x ) 


= C-4 


—1 

1 


C-i can t> e constrained by equating xo and Co a t x = — (j due to the continuity of the wavefunction: C_i = —XS) X-'- = *Xs 
where xo is a real constant to be fixed by the normalization condition. 

The 4-component wavefunction of the left Majorana zero mode T F thus becomes (up to normalization), 

( (expfi£(a;+ §)) ) 

( * 


&l(x) = XSs < 


, ( i v^( X +§ 

- ex p(-*£( x +i)) 

~ ex p(^( x + §)) 

[ie x p(-i£(x + §)) ) 


e ^ x+ 2 ^0(x + —) + 


-1 

-1 

i 


e v F ( x +2 )q(_ x _ —) . (B5) 
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Right Majorana zero mode. Since 'H[J,x) = K'Hr{—x)K, we can directly derive the wavefunction for the right Majorana 
zero mode 'l'fl(x) = K^l(-x) = x). 

Zero-bias splitting. If we assume |<3>+(x)) is the eigenstate of the full BdG Hamiltonian 'H(x) with eigenenergy £ + , 
namely 7~L (x)| < f>+(x)) = £+|<l>+(ir)), then according to the particle-hole symmetry, we will have 7~L(x) (S|<J>+(x))) = 
—£ + (S|<l>+(x))), therefore we denote |(cc)) = 5|$+(x)), which satisfies 'H(x)|<I , “(x)) = £_|<&“(x)} = 

—£+\$~g (x)). One can check that S|\I/^(x)) = ^/{(x)), S| v I'l(x)) = | v ki(x)), thus we can safely assume (a:)) = 
^ (I^Mx)) ± t|*i(x))), |$±(-X)) = ± (l^fl(-x)) ± i|$ L (-x))) = ±%K |<b±(x)>. 

The eigenenergies of these two approximate eigenstates are given by £± « =p Im / 0 °° ^' r (x)1Ar{x)^! i(x)dx , which can be 
further simplified: / 0 +?c '&] i (x)Un(x)'I/ l(x)cIx = Wf^r{x = 0)(t 2 <8> <t z )T'l(x = 0), where we have integrated the kinetic 
term by parts. Therefore, we derive: £± sa =p x/.-Rc 


^ R {x = 0)(r 2 ® ct 2 )'T l (x = 0) 


Thus far, our analysis is completely general, having used only the symmetries of the problem at hand. To obtain the final 
expression, we plug in the explicit forms for ^l/r, which gives: 


£ ± = ^ 


-AL / fiL 


± 2 vfXq e vp cos ( — ) = i 


A B _m / /rL 


A + B 


’“F COS 


V F 


B» A -fi f /rL\ 

->• ± A e v f cos — . 

\vf ) 


(B 6 ) 


A parallel evaluation can be worked out for the SC-FM-SC setup with the inclusion of a finite magnetic field term. 

It is instructive to compare this result with that obtained in the main text. Importantly, we find that the oscillatory cosine term 
has exactly the same form in both calculations, and vanishes if the chemical potential is tuned to zero. However, they differ in 
the exponential decay rate. Let us summarize the key result of this section: the zero-energy splitting in the quantum spin-Hall 
edge has an oscillatory dependence on //, which exactly matches that found in the instanton calculation. 


Appendix C: Instantons in finite-interval sine-Gordon model with Neumann boundary conditions 


In the instanton calculation presented in the main text, a crucial element is the evaluation of the Fredholm determinant. For 
periodic boundary conditions, this has been discussed for the sine-Gordon theory in Ref. 49 ; here we will review the relevant 
aspects of the derivation, and explain how the situation differs for open boundaries. 

To evaluate the functional determinant, we must first explicitly factor out the zero-mode resulting from time-translation of 
the center of the instanton. This leads to a factor of So/(2nH) (see, e.g. Ref. 55 ). In addition, one typically evaluates a 
ratio of the functional determinant of interest to the functional determinant of a harmonic oscillator, including the latter into the 
normalization factor. 

Using the dilute instanton gas approximation, in this approach formally the energy splitting is given by 49 




—2hy/vcosf—q<f, - 

\n hnvJ 





det + t 


^/det' [-dl, - <9 2 , - ^ sin( 2 ru/w(T))] 


(Cl) 


where we have rescaled the variables and the fields, and the prime on det 7 means that the zero eigenvalue has been excluded. 

In a one-dimensional system, the ratio of determinants in Eq. (Cl) is divergent due to spatial fluctuations, and must be 
regulated. As shown in Ref. 56 , the divergent contribution to the functional determinant ratio is canceled by using renormalized 
values of the effective mass and coupling constant in So, and the finite part can be expressed: 


, f det' M \ d , d ( 1 f 0 

log! TTXv 1 — = -- (p^ 


dt t z ~ 1 {K t (M) - K t (M 0 ) - 1} 


v detM 0 y dz" ’ dz \L\z) Jo ' " ' " '/ *->-o 

Here the heat kernel K t (A) = Tr[e _t " 4 ] and the Hermitian operators M, M 0 and (). Qq are defined as follows. 


M = -d 2 x , - d 2 T , + ml ( 1 - 


cosh 2 (to 0 (t' - t£)) 


= - d 2 , + Q, M 0 = -d 2 x , - dl, + TOq = —dl, + Q 0 , 


(C2) 


(C3) 


where the mass term mo = \/4 A/(Ir). 

It is useful to separate the zeta function into two main pieces: £( 2 , M) = (i 0 g(z, M) + Cmix(~j M), where 
Clog (Z, M) = ^ J™ dt t z ~ l { I< t {-dl,) - 1} , Cmix(*, M) = Y^j^ dt tZ ~ l { K t(-d 2 x .) ■ [l<t(Q) - K t (Qo) - l] } , 
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with Q , Q o given in Eq. (C3). _ 

Of primary interest is the contribution of £i og , which gives rise to the extra l/y/L (L = yfvL 1 ) factor in the prefactor of the 
energy splitting, leading to an T-independent prefactor in agreement with known results for Majorana zero modes 33,35 . This 
term can be computed exactly: 


Ciog( 2: ) -^0 — 


r 0 ) 


dff z - 1 {iT t (-e)-l} = E 




dt t z 


In the continuum limit with the von Neumann boundary conditions, the momentum k n = 



Clog(~; M) 



( L' 2 \ z 

f^-J c(2z) and exp 


1 ^Ciog A 

2 dz z= o J 


1 1 _ vi 1 


Crucially, here we have kept only positive k values in the sum, as is appropriate for general (i.e. fixed or free) open boundary 
conditions. For periodic boundary conditions, where both positive and negative k are allowed, Clog is proportional to 1/ L 56 . 

The remaining term Cmix is important chiefly because, for open boundary conditions, it contributes a dimensionful constant to 
the prefactor. It is convenient to divide this into the following three contributions: 


Cmix(z, M) 
Ccon (Z, M) 


Clin (Z,M) 


Cexp {Z, M) 


CconZ, M) + Clin M) + Cexp(~, M), 


1 


T(z) 

1 

m 

1 


dt t z 


1 


K t {Q) - K t (Q 0 ) - 1 


1 ( T' 

( 4V47T t 


r °° dtt^UKti-dl,)- ^ 


K t (Q) - K t (Qo) ^ 1 
2 4-^Ant J 


K t (Q) - K t (Q 0 ) 



Cexp (z,M) gives an additive contribution that is exponentially suppressed in L relative to the other two terms, and can be 
neglected. We will show that the Ccon-contribution - which is not present with periodic boundary conditions - ensures the correct 
dimensionality for the splitting. 

To compute these, we require the spectra and densities of states of the operators Q and Qo- Apart from the zero-mode of Q, 
which we have already explicitly factored out, these are continuous functions, and are given in Ref. 49 . This gives: 


Ccon(2, M) 


1 

m 



1 

2 


- r d P ( 2m ° 2 ) 

nJ 0 \P 2 + m oJ 


e ~t(p 2 +ml) 


1 - 2 z r(z + i/ 2 ) 

2 \dn 0 T(z +1 ) 


Therefore, the extra contribution to the prefactor assumes 


1 d(con(z,M) 


2 dz 


z =0 


1 

2 


log (2mo) and exp 


1 d( con (z, M) 

2 dz 





(C4) 


where we use the definition of ui = 2yJ A v/(H<;) = Too\/^ i n the main text. This ensures that the overall prefactor has the 
dimensions of energy. For periodic boundary conditions, where the contribution of Clog ~ 1/L, there is no Ccon-contribution. 

Ciin can be calculated in a similar manner, and makes an additive contribution to So that scales linearly with L. This must be 
combined with other such finite terms that arise in fixing the renormalized value of m relative to the bare value appearing in Sq. 
(Recall that it is necessary to use this renormalized value, as it cancels the divergences in the determinant ratio in Eq. (Cl) 56 ). 
Ref. 49 verified numerically that the net effect of these finite contributions on the energy splitting is negligible for the sine-Gordon 
model with periodic boundary conditions; since these contributions are independent of the boundary conditions, we expect that 
this holds in the present context as well, and leave a more precise evaluation of the leading quantum corrections to Sq for the 
future investigations. 

The final finite result for the energy splitting thus assumes the following form. 




2 tiui 
n 



7T /iL 

n ^ Hnv 


(C5) 


After incorporating both the spatial and temporal fluctuations, the overall dependence of the prefactor on the system’s size 
becomes independent of L, which results from a cancelation between the factor of t/Sq due to the time-translational-invariant 
zero-mode, and the effect of finite-momentum fluctuations in the Fredholm determinant with the Neumann boundary conditions. 

































10 


* chen2698@umn.edu 
t fburnell@umn.edu 

1 A. Y. Kitaev, Phys.-Usp. 44, 131 (2001). 

2 L. Fu and C. L. Kane, Phys. Rev. B 79, 161408 (2009). 

1 R. M. Lutchyn, J. D. Sau, and S. Das Sarma, Phys. Rev. Lett. 105, 
077001 (2010). 

4 Y. Oreg, G. Refael, and F. von Oppen, Phys. Rev. Lett. 105, 
177002 (2010). 

5 S. R. Elliott and M. Franz, Rev. Mod. Phys. 87, 137 (2015). 

6 P. Fendley, J. Stat. Mech. 2012, 11020 (2012). 

7 M. Barkeshli and X.-L. Qi, Phys. Rev. X 2, 031013 (2012). 

8 D. J. Clarke, J. Alicea, and K. Shtengel, Nat. Commun. 4, 1348 
(2013). 

9 N. H. Lindner, E. Berg, G. Refael, and A. Stern, Phys. Rev. X 2, 
041002 (2012). 

10 M. Cheng, Phys. Rev. B 86, 195126 (2012). 

11 A. Vaezi, Phys. Rev. B 87, 035132 (2013). 

12 R. S. K. Mong, D. J. Clarke, J. Alicea, N. H. Lindner, P. Fendley, 
C. Nayak, Y. Oreg, A. Stern, E. Berg, K. Shtengel, and M. P. A. 
Fisher, Phys. Rev. X 4, 011036 (2014). 

13 D. J. Clarke, J. Alicea, and K. Shtengel, Nat. Phys. 10, 877 (2014). 

14 J. Klinovaja and D. Loss, Phys. Rev. Lett. 112, 246403 (2014). 

15 A. S. Jermyn, R. S. K. Mong, J. Alicea, and P. Fendley, Phys. Rev. 
B 90, 165106(2014). 

16 M. Barkeshli, C.-M. Jian, and X.-L. Qi, Phys. Rev. B 88, 235103 
(2013). 

17 V. Mourik, K. Zuo, S. M. Frolov, S. R. Plissard, E. P. A. M. 
Bakkers, and L. P. Kouwenhoven, Science 336, 1003 (2012). 

18 A. Das, Y. Ronen, Y. Most, Y. Oreg, M. Heiblum, and H. Shtrik- 
man, Nat. Phys. 8, 887 (2012). 

19 M. T. Deng, C. L. Yu, G. Y. Huang, M. Larsson, P. Caroff, and H. 
Q. Xu, Nano Lett. 12, 6414 (2012). 

20 H. O. H. Churchill, V. Fatemi, K. Grove-Rasmussen, M. T. Deng, 
P. Caroff, H. Q. Xu, and C. M. Marcus, Phys. Rev. B 87, 241401 
(2013). 

21 S. Hart, H. Ren, T. Wagner, P. Leubner, M. Miihlbauer, C. Briine, 
H. Buhmann, L. W. Molenkamp, and A. Yacoby, Nat. Phys. 10. 
638 (2014). 

22 V. S. Pribiag, A. J. A. Beukman, F. Qu, M. C. Cassidy, C. Charp- 
entier, W. Wegscheider, and L. P. Kouwenhoven, Nat. Nanotech- 
nol. 10, 593 (2015). 

2j A. P. Higginbotham, S. M. Albrecht, G. Kirsanskas, W. Chang, F. 
Kuemmeth, P. Krogstrup, T. S. Jespersen, J. Nygard, K. Flensberg, 
and C. M. Marcus, Nat. Phys. 11, 1017 (2015). 

24 A. D. K. Finck, C. Kurter, E. D. Huemiller, Y. S. Hor, D. J. Van 
Harlingen, arXiv:1503.06898vl (2015). 

25 S. Nadj-Perge, I. K. Drozdov, J. Li, H. Chen, S. Jeon, J. Seo, A. 
H. MacDonald, B. A. Bernevig, and A. Yazdani, Science 346, 602 
(2014). 

26 J. Alicea, Y. Oreg, G. Refael, F. von Oppen, and M. P. A. Fisher, 
Nat. Phys. 7, 412(2011). 

27 M. H. Freedman, A. Kitaev, M. J. Larsen, and Z. Wang, Bull. 
Amer. Math. Soc. 40, 31 (2003). 

28 A. Y. Kitaev, Ann. Phys. (N.Y.) 303, 2 (2003). 

29 C. Nayak, S. H. Simons, A. Stern, M. Freedman, and S. Das 
Sarma, Rev. Mod. Phys. 80. 1083 (2008). 

30 S. Das Sarma, M. Freedman, and C. Nayak, Phys. Rev. Lett. 94, 
166802 (2005). 

31 P. Bonderson, Phys. Rev. B 87, 035113 (2013). 


32 J. Alicea and P. Fendley, arXiv:1504.02476vl (2015). 

33 S. Das Sarma, J. D. Sau, and T. D. Stanescu, Phys. Rev. B 86, 
220506 (2012). 

34 L. Fidkowski. R. M. Lutchyn, C. Nayak, and M. P. A. Fisher, 
Phys. Rev. B 84, 195436 (2011). 

35 A. A. Zyuzin, D. Rainis, J. Klinovaja, and D. Loss, Phys. Rev. 
Lett. Ill, 056802 (2013). 

36 D. Rainis, L. Trifunovic, J. Klinovaja, and D. Loss, Phys. Rev. B 
87, 024515 (2013). 

’ 7 F. Maier, J. Klinovaja, and D. Loss, Phys. Rev. B 90, 195421 
(2014). 

38 M. Cheng, R. M. Lutchyn, V. Galitski, and S. Das Sarma, Phys. 
Rev. Lett. 103, 107001 (2009). 

39 M. Cheng, R. M. Lutchyn, V. Galitski, and S. Das Sarma, Phys. 
Rev. B 82, 094504 (2010). 

40 M. Baraban, G. Zikos, N. Bonesteel, and S. H. Simon, Phys. Rev. 
Lett. 103, 076801 (2009). 

41 V. Lahtinen, G. Kells, A. Carollo, T. Stitt, J. Vala, and J. K. Pa- 
chos, Ann. Phys. (N.Y.) 323, 2286 (2008). 

42 R. S. K. Mong, D. J. Clarke, J. Alicea, N. H. Lindner, and P. Fend¬ 
ley, J. Phys. A: Math. Theor. 47 , 452001 (2014). 

43 W. Li, S. Yang, H.-H. Tu, and M. Cheng, Phys. Rev. B 91, 115133 
(2015). 

44 E. M. Stoudenmire, D. J. Clarke, R. S. K. Mong, and J. Alicea, 
Phys. Rev. B 91, 235112 (2015). 

45 Y. Zhuang, H. J. Changlani, N. M. Tubman, and T. L. Hughes, 
Phys. Rev. B 92, 035154 (2015). 

46 X.-G. Wen, Quantum Field Theory of Many-Body Systems, (Ox¬ 
ford University Press, New York, 2004). 

47 X.-G. Wen, Advances in Physics 44, 405 (1995). 

48 T. Giamarchi, Quantum Physics in One Dimension, (Oxford Uni¬ 
versity Press, New York, 2004). 

49 Z. Bajnok, L. Palla, G. Takacs, and F. Wagner, Nucl. Phys. B 587, 
585 (2000). 

50 Y. I. Frenkel and T. K. Kontorova, Zh. Eksp. Teor. Fiz. 8, 1340 
(1938) [Phys. Z. Sowjetunion 13, 1 (1938)]. 

51 T. H. R. Skyline, Proc. Roy. Soc. Lond. A 247 , 260 (1958). 

52 T. H. R. Skyrme, Proc. Roy. Soc. Lond. A 262 , 237 (1961). 

53 A. M. Polyakov, Nucl. Phys. B 120, 429 (1977). 

54 S. Coleman, Aspects of symmetry, (Cambridge University Press, 
New York, 1985). 

55 A. I. Vainshtein, V. I. Zakharov, V. A. Novikov, and M. A. Shif- 
man, Usp. Fiz. Nauk 136 , 553 (1982) [Sov. Phys. Usp. 25, 195 
(1982)]. 

56 G. Munster, Nucl. Phys. B 324 , 630 (1989). 

57 V. Shivamoggi, G. Refael, and J. E. Moore, Phys. Rev. B 82, 
041405 (2010). 

58 P. Bonderson, Phys. Rev. Lett. 103, 110403 (2009). 

59 M. Barkeshli, P. Bonderson, M. Cheng, and Z. Wang, 
arXiv:1410.4540v2 (2014). 

60 P. Fendley, J. Phys. A: Math. Theor. 47, 075001 (2014). 

61 R. J. Baxter, J. Phys. A: Math. Theor. 47, 315001 (2014). 

62 H. Au-Yang and J. H. H. Perk, J. Phys. A: Math. Theor. 47, 
315002 (2014). 

63 It is worth noticing that the chiral phases are not allowed in the 
Ising chain. 

64 E. Fradkin and L. P. Kadanoff, Nucl. Phys. B 170. 1 (1980). 



